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1 Introduction 



1.1 Bayes linear methods 

Bayes linear methods are a form of Bayesian statistics, which acknowledge the difficulties asso- 
ciated with the full modelling, specification, and conditioning required by distributional Bayesian 
statistics, and instead try to make best possible use of partial specifications, based on means, vari- 
ances and co variances. Unsurprisingly, much of the theory is formally identical to inference in 
multivariate Gaussian Bayesian networks, but interpretation of results is generally different. This 
document assumes a working knowledge of the basic tools of the Bayes linear methodology. An 
introduction to Bayes linear methods is given in [[]]]. An introduction to (non-local) computational 
issues can be found in [§]. The foundations of the theory are discussed in [Qj, [|3|], and [§]. On-line, 
an introduction to the theory can be found in [f|], from the Bayes Linear Methods WWW home 
page: http://fourier.dur.ac.uk:8000/ stat s /bayes 1 in/ 

1.2 LISP-STAT 

LISP-STAT is an interpreted, object-oriented environment for statistical computing, described in 
|Q] . This document assumes a working knowledge of LISP-STAT, and the basics of object-oriented 
programming. On-line, LISP-STAT information is available from the LISP-STAT WWW home 

page: http : / / www . stat . umn . edu/ ~ luke/ xls/xlsinfo/xlsinfo. html 

1.3 Local computation 

BAYES-LIN carries out local computation via message-passing between adjacent nodes of a clique- 
tree representing the statistical model of interest. Again, local computation in Bayesian networks 
is a huge area, and this document assumes a working knowledge of graphical models, conditional 
independence and some of the ideas behind local computation. The best introduction to all of these 
areas is [0] . In particular, Chapter 3 of that volume deals with all of the relevant graph-theoretic 
concepts, and Section 7.2 gives an introduction to graphical Gaussian models. 

1.4 Installing and running BAYES-LIN 

You need a working LISP-STAT system installed before you attempt to install BAYES-LIN. The 
following instructions are for a UNIX system with an XLISP-STAT installation, but installing on 
other systems should be similar. Note that the graphics work best on systems with at least a 16 bit 
colour display. If you only have an 8 bit display (256 colours), make sure that most are free for 
use by BAYES-LIN. The graphics will not work on displays poorer than 8 bit colour. Create a new 
directory for the BAYES-LIN system. Download the BAYES-LIN software from the BAYES- 
LIN WWW page: http : / / www . ncl . ac . uk/ ~nd jwl /bayes 1 in/ and put into the new 
directory. In this new directory type: 

% gunzip blinO la . tar . gz 
% tar -xvf blinOla.tar 
% gzip blinOla.tar 

You should then be able to run LISP-STAT with the BAYES-LIN extensions simply by running 

% xlispstat 
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from within this directory. You can check that the extensions are loaded by typing in some of the 
following commands in the LISP-STAT listener window. 

> (help ' create-tree-node) 

> (send moral-node-proto :help) 

> (send tree-node-proto :help : observe) 

In general, to make sure the extensions are loaded, use the expression 

> (require "bayeslin") 

When you are satisfied that the extensions are loaded, exit BAYES-LIN. 

> (exit) 

In order to run the examples, simply call LISP-STAT with the example as first argument, eg. 

% xlispstat ex-dlm 
or 

% xlispstat ex-mdlm 

These two examples will be explained in the following sections. 



2 A "toy" dynamic linear model 



2.1 Description of the model 

The BAYES-LIN code for this example can be found in the file ex-dlm . 1 sp, which is part of the 
standard BAYES-LIN distribution. The example concerns a very simple model for 3 observations 
in time. The model can be written in the form of a locally constant DLM. 

X t = 9 f +V f 

e r = e r _i+co f 

X t denotes the observation at time t (t = 1,2, 3), which is dependent on the state of the system at 
time t, Q t . The variables v t and (O t are incidental noise terms. The model is initialised by specifying 
beliefs about the initial state of the system; in this case, E (0i) = 1, Var (9i) = 1, and the variance 
of the noise terms; in this case, Var (co f ) = Var (v f ) = 1 . Of course, for such a simple model, a non- 
local analysis is trivial, since the expectation vector for the entire system, and the variance matrix 
for the entire system can be written down and worked with directly. E (0i, 02, 03,Xi,X2,X3) = 
(1,1,1,1,1,1) T , 

/ 1 1 1 1 1 1 \ 



Var(0 1 ,0 2 ,03,Z 1 ,X 2 ,Z 3 ) = 



V i 



1 2 2 

1 2 3 

1 1 1 

1 2 2 

2 3 



1 
1 
2 
1 
1 



2 2 

2 3 

1 1 

3 2 

2 4/ 



However, for the example in the next section, such explicit non-local analysis will not be possible. 



4 



2.2 Graphical models 

The graph for this model is shown below. 




Since there are no unmarried parents, this graph can be moralised simply by dropping arrows. 



mT1 rnT2 rnT3 




mX1 mX2 mX3 



Now, since there are no cycles, this graph is already triangulated, so the clique tree may be 
formed as follows. 

tT12 tT23 




tX1 tX2 tX3 



2.3 Defining the clique-tree 

BAYES-LIN carries out computation on the clique-tree, and displays results on the moral graph. 
Therefore, both need to be introduced to the BAYES-LIN system. Since all information and com- 
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putations are carried out on the clique-tree (in fact, computations can be carried out without defin- 
ing a moral graph at all), this is defined first. Appropriate code for defining the tree nodes is shown 
below. 

(create-tree-node ' tXl ' (tl xl) #(1 1) #2a((l 1) (1 2)) ' (tT12) ) 
(create-tree-node ' tX2 ' (t2 x2) #(1 1) #2a((2 2) (2 3)) ' (tT23) ) 
(create-tree-node ' tX3 ' (t3 x3) #(1 1) #2a((3 3) (3 4)) ' (tT23) ) 
(create-tree-node 'tT12 ' (tl t2) #(1 1) #2a((l 1) (1 2)) ' (tXl tT23)) 
(create-tree-node 'tT23 ' (t2 t3) #(1 1) #2a((2 2) (2 3)) ' (tT12 tX2 tX3)) 

The global function create-tree-node is used to define each node in turn. The function 
expects five arguments. The first argument is a symbol to point to the resulting tree-node object. 
The second is a list of variables which the node contains. The third and fourth are the expectation 
vector and variance matrix for the variable list, and the fifth is a list of neighbouring tree nodes. 
Next, the moral nodes are defined. 



2.4 Defining the moral graph 



(create-moral-node 'mXl '(xl) ' tXl "mXl" '(mTl)) 

(create-moral-node ' mX2 '(x2) ' tX2 "mX2" ' (mT2) ) 

(create-moral-node 'mX3 ' (x3) ' tX3 "mX3" ' (mT3) ) 

(create-moral-node 'mTl '(tl) ' tXl "mTl" ' (mXl mT2)) 

(create-moral-node 'mT2 ' (t2) ' tX2 "mT2" ' (mX2 mTl mT3) 

(create-moral-node 'mT3 ' (t3) ' tX3 "mT3" ' (mX3 mT2)) 



The global function create-moral-node is used to define each node in turn. The first is 
a symbol to bind the object to. The second is a variable list. The third is a clique-tree node which 
contains all of the variables at this node (such a node always exists). The fourth is a string to be 
used for plotting purposes, and the fifth is a list of neighbouring moral graph nodes. Next, some 
plotting positions are defined by sending a :location message to each moral node object. 

(send mXl rlocation '(0.2 0.8)) 
(send mX2 :location '(0.5 0.8)) 
(send mX3 :location '(0.8 0.8)) 
(send mTl :location '(0.2 0.2)) 
(send mT2 rlocation '(0.5 0.2)) 
(send mT3 :location '(0.8 0.2)) 

This step may be omitted if one is not interested in plotting of results. The locations are on a (0, 1) 
scale for x and y coordinates, respectively. The origin is the top-left of the plot window. The model 
is now completely specified. Before carrying out adjustment, we create plot windows to show 
diagnostic information. 



2.5 Adjustments 

(create-moral-plot 'myplot) 

(ere ate -global -moral -plot ' myplot 2 ) 
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This creates a plot window with the name myplot to show partial adjustment information, and 
another, myplot 2, to show global adjustment information. Note that high-quality colour Encap- 
sulated PostScript output is produced for each plot after each redraw of the screen, and stored in 
the files mpw . eps and gmpw . eps respectively. 

We are now in a position to carry out adjustments. Suppose that variable X\ is observed to be 
x\ = 3. This information can be introduced into the graph by sending the following message to the 
appropriate moral graph node. 

(send mxl robserve ' (xl) #(3)) 

In general, one can observe a list of variables, provided all variables are contained in the moral 
graph node receiving the message. The message is passed on to the appropriate clique-tree node, 
and then propagated around the clique-tree. We can tell our plot object to gather information from 
the tree for display, as follows. 

(send myplot : record) 

The plot should now show how information flows around the moral graph (more on this later). Note 
that although information has been introduced into the graph, it has not been absorbed into it, and 
that further information can not be introduced until it has. This can be understood by sending some 
messages to the graph, and looking at the return values. If the expectation and variance of the first 
observable node is examined 

(send mxl rex) 
(send mxl rvar) 

it can be seen that it retains its a priori values. However, one can also ask for adjusted expectations 
and variances. 

(send mxl raex) 
(send mxl ravar) 

Similar queries can be sent to the third moral graph node. 

(send mx3 rex) 
(send mx3 rvar) 
(send mx3 raex) 
(send mx3 ravar) 

When we are finished examining the effects of the current adjustment, and wish to add further 
information into the graph, the current information should be absorbed. 

(send mxl r absorb) 

The absorbing makes the adjusted information the new prior information, ready for the next adjust- 
ment. This can be verified by looking at the new expectation and variance for the third observation. 

(send mx3 rex) 
(send mx3 rvar) 

Finally, we can introduce new information, record it, and then absorb it, before examining the 
results. 
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(send mx2 :observe ' (x2) #(-1)) 
(send myplot : record) 
(send mx2 : absorb) 

(send mx3 :var) 
(send mx3 rex) 

The two plot windows should now look similar to the following. 
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2.6 Interpreting the graphics window 

Whenever an observation is made and recorded, a portion is removed from the outside of each 
node. The area removed is proportional to the variance resolved by the adjustment. Consequently, 
the radius removed is proportional to the standard deviation resolved. Therefore, when a node is 
fully observed, there is no dark centre left remaining. For other nodes, the size of the dark centre 
is proportional to the proportion of original uncertainty left remaining. For multivariate nodes, the 
Bayes linear concept of resolution is used. 

The additional red and blue shadings give an indication of the changes in expectation, relative 
to a priori uncertainty. Red shadings indicate changes larger than expected, and blue shadings 
represent changes in expectation smaller than expected^. The amount of red and blue shading 
increases as the "degree of surprisingness" increases. The amount of shading is a transformation 
of the Bayes linear concept of size-ratio. The transformation can be user-specified by redefining 
the plot-object's : sr-map method appropriately. 



3 Computation for a large multivariate DLM 
3.1 Description of the model 

The BAYES-LIN code for this example can be found in the file ex-mdlm . lsp, which is part of 
the standard BAYES-LIN distribution. The following data represent weekly sales of six soft drinks 
packs from a wholesale depot. 

51 27 1 4 6 3 
113 55 7 15 4 
103 71 10 16 7 



Clearly a multivariate time series model is required for such data. The following multivariate 
locally constant DLM is adopted. 

X t = Q t +V t 
9/ = 9 f -i +a>f 

This is the same model as used in the last example, but here all of the variables denote random 
vectors of dimension six. The model is specified in the following way. There are 35 observa- 
tions, and so t runs from 1 to 35 for the actual observations. However, for this model, it was felt 
more convenient to initialise the model at t = 0. The initial state of the system was specified as 
E (6 ) = (50, 50, 50, 50, 50, 50) T and Var (0 O ) = diag(900, 900, 900, 900, 900, 900) . The covariance 

*Note that changes in expectation smaller than expected can still be of concern, since they are indicative of a 
possible under-utilisation of prior information. 
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structure for the noise terms was specified to be Var (v f ) = V, Var ((o f ) = W, where 



V = 



W = 



/ 
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2420 
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4.78 


1 

\ 


Jo 


D 1 


j 


97 


c 

J 


QQ 


1 1 ^ 

i i.jj 


4 78 
4-. / o 


44.21 


J 


/ 


1112 


49 


272 


47 


22 


52 


66.45 


31.56 


27.84 \ 






272 


47 


195 


50 


11 


53 


30.07 


18.51 
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07 
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54 
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14.04 


8.03 






31 


56 


18 


51 


4 
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14.04 


40.50 


7.32 




V 


27 


84 


15 


37 


6 


28 


8.03 


7.32 


32.97 J 





See for an explanation of the given specification. These specifications determine the model, but 
note that there are6x4x35 + 6 = 846 variables in this problem (assuming that we are interested 
in the noise terms). This problem is about at the limit of the size which can be tackled by a brute 
force approach, making a local computation approach particularly attractive. 

3.2 Graphical models 

We are interested in making inferences about the noise terms in this example (in order to help 
diagnose deficiencies of the model), and so the noise terms need to be included in the model. The 
first part of the DAG for this structure is therefore as follows (note that the DAG nodes are all 
multivariate). 




Marrying parents and dropping arrows gives the moral graph for the problem (note that the 
moral graph nodes are all multivariate). 
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3.3 Some constants 

The code for such a problem can be constructed as follows, 
matrix, and some constants are defined. 



First, the data is read, turned into a 



(def mydata (read-data-columns "ex-mdlm . dat " ) ) 

(def data (make-array '(6 35) : initial-contents mydata)) 



(def v 
#2a( 



2420.36 
387 . 330 
20.3907 
165.274 
44.5645 
58 . 6081 



387.330 
263.850 
3 . 85480 
71 . 5054 
23.4794 
3.26543 



20.3907 
3 . 85480 
30 .7910 
3 . 58376 
1.25836 
5.98943 



165.274 
71 . 5054 
3 . 58376 
139.715 
23 . 1193 
11 . 3268 



44.5645 
23.4794 
1.25836 
23 . 1193 
50 . 0087 
4 .78345 



58 . 6081 
3 .26543 
5 . 98943 
11 . 3268 
4 .78345 
44 . 2135 



(def w 
#2a ( 

( 1112.49 

( 272.473 

( 22.5176 

( 66.4472 

( 31.5611 

( 27.8440 



272 .473 
195.499 
11.5298 
30.0701 
18 . 5065 
15 . 3737 



22 . 5176 
11.5298 
29 . 6411 
5.53985 
4 .66797 
6.27591 



66.4472 
30.0701 
5.53985 
78.9076 
14 . 0421 
8.03411 



31.5611 
18.5065 
4 .66797 
14 . 0421 
40 . 5035 
7.32038 



27 . 8440 
15.3737 
6.27591 
8 . 03411 
7.32038 
32 . 9728 
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(def eO (coerce (repeat 50 6) 'array)) 

(def eeO (coerce (append (repeat 50 12) (repeat 6)) 'array)) 
(def wO (diagonal (repeat 900 6) ) ) 
(def zero66 (diagonal (repeat 6) ) ) 

These are all self-explanatory. 

3.4 Defining the clique-tree 

We can now create the type B cliques as follows. 



(dolist (i (iseq 1 35)) 
(create-tree-node 
(intern (format nil "b~a" i) ) 



(list 


( intern 


(format nil 


"xl 


~a" 


( intern 


(format 


nil 


"x2 . ~a" 


i) ) 




( intern 


(format 


nil 


"x3 . ~a" 


i) ) 




(intern 


(format 


nil 


"x4 . "a" 


i) ) 




(intern 


(format 


nil 


"x5 . "a" 


i) ) 




( intern 


( format 


nil 


" x 6 . " a " 


i) ) 




( intern 


( format 


nil 


"thetal . 


"a" 


i) ) 


( intern 


( format 


nil 


"theta2 . 


"a" 


i) ) 


( intern 


( format 


nil 


"theta3 . 


"a" 


i) ) 


( intern 


( format 


nil 


"theta4 . 


"a" 


i) ) 


( intern 


( format 


nil 


"theta5 . 


"a" 


i) ) 


( intern 


( format 


nil 


"theta6 . 


"a" 


i) ) 


( intern 


( format 


nil 


"nul . ~a' 


i ) 




( intern 


( format 


nil 


"nu2 . ~a' 


i ) 




(intern 


(format 


nil 


"nu3 . ~a' 


i ) 




( intern 


(format 


nil 


"nu4 . ~a' 


i ) 




( intern 


(format 


nil 


"nu5 . ~a' 


i ) 




(intern 


(format 


nil 


"nu6 . "a' 


i ) 





eeO 

(bind-rows (bind-columns 
(+ wO (* w i) v) 
(+ wO (* w i) ) 
v) 

(bind-columns 

(+ wO (* w i) ) 

(+ wO (* w i) ) 

zero66) 
(bind-columns 

v 

zero66 
v) 

) 

(if (= i 35) 

(list (intern (format nil "a~a" i) ) ) 

(list (intern (format nil "a~a" i) ) 
(intern (format nil "a~a" (+ i 1) ) ) ) 

) 

) ) 

Note that the expression (intern (format nil "b~a" i)) means "create the Lisp 
symbol hi, where i is a variable". This trick is used a lot for the construction of big models with a 
repetitive structure. Next, the type A cliques can be constructed, in a very similar way. 



(dolist (i (iseq 1 35)) 
(create-tree-node 



(intern (format nil "a~a" 


i) ) 








(list (intern 


( format 


nil 


"thetal . 


"a" ( 


- i 


1) ) 


( intern 


( format 


nil 


"theta2 . 


~a" ( 


- i 


1) ) 


( intern 


( format 


nil 


"theta3 . 


~a" ( 


- i 


1) ) 


( intern 


( format 


nil 


"theta4 . 


~a" ( 


- i 


1) ) 


( intern 


( format 


nil 


"theta5 . 


~a" ( 


- i 


1) ) 


( intern 


( format 


nil 


"theta6 . 


~a" ( 


- i 


1) ) 


( intern 


( format 


nil 


"thetal . 


"a" i 


) ) 
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(intern 


(format nil 


"theta2 


"a" 


i) ) 




(intern 


(format nil 


"theta3 


~a" 


i) ) 




(intern 


(format nil 


"theta4 . 


~a" 


i) ) 




(intern 


(format nil 


"theta5 . 


"a" 


i) ) 




(intern 


(format nil 


"theta6 . 


"a" 


i) ) 


(intern 


( format 


nil "omegal 


~a" i) ) 






( intern 


( format 


nil "omega2 


~a" i) ) 






( intern 


( format 


nil "omega3 


~a" i)) 






(intern 


( format 


nil "omega4 


~a" i) ) 






(intern 


( format 


nil "omega5 


~a" i) ) 






( intern 


( format 


nil "omega6 


~a" i) ) 







) 

eeO 



(bind-rows (bind-columns 
(+ wO (* w (- i 1) ) ) 
(+ wO (* w (- i 1) ) ) 
zero66) 
(bind-columns 
(+ wO (* w (- i 1) ) ) 
(+ wO (* w i) ) 
w) 

(bind-columns 
zero66 
w 
w) 

) 

(if (= i 1) 

(list (intern (format nil "b~a" i) ) ) 
(list (intern (format nil "b~a" i) ) 
(intern (format nil "b"a" (- i 1)))) 
) 

) ) 



3.5 Defining the moral graph 

Next, moral graph nodes need to be created, for diagnostic plotting purposes. Since there isn't 
room on the average computer screen for the moral graph for all 35 time point, the structure will 
only be constructed for the first 8 time points only. 

; ; number of moral nodes to create and plot 
(def plotnum 8) 

(create-moral-node (intern "theta.O") 

(list (intern "thetal.O") 
(intern "theta2.0") 
(intern "theta3.0") 
(intern "theta4.0") 
(intern "theta5.0") 
(intern "theta6.0")) 

(intern "al") 

"Theta (0) " 

(list (intern "theta. 1") 
(intern " omega. 1")) 
) 

(send (symbol-value (intern "theta.O")) :location (list (/ 1 (+ plotnum 2)) 0.4)) 
(dolist (i (iseq 1 plotnum) ) 
; ; create the theta node 

(create-moral-node (intern (format nil "theta. "a" i)) 

(list (intern (format nil "thetal . "a" i) ) 
(intern (format nil "theta2."a" i)) 
(intern (format nil "theta3.~a" i)) 
(intern (format nil "theta4 . ~a" i)) 
(intern (format nil "theta5.~a" i)) 
(intern (format nil "theta6.~a" i) ) ) 
(intern (format nil "b~a" i)) 
(format nil "Theta (~a)" i) 
(if (< i plotnum) 
(list (intern (format nil "omega. ~a" i) ) 
(intern (format nil "nu.~a" i)) 
(intern (format nil "x.~a" i)) 
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(intern (format nil "theta.~a" (- i 1))) 
(intern (format nil "theta.~a" (+ i 1))) 
(intern (format nil "omega. "a" (+ i 1))) 
) 

(list (intern (format nil "omega. "a" i) ) 
(intern (format nil "nu. "a" i)) 
(intern (format nil "x.~a" i)) 
(intern (format nil "theta.~a" (- i 1))) 
) 

) 

) 

(send (symbol-value (intern (format nil "theta.~a" i) ) ) :location 
(list (* (+ i 1) (/ 1 (+ plotnum 2))) 0.4) 

) 

; ; create the omega node 

( create-moral-node (intern (format nil "omega. "a" i) ) 

(list (intern (format nil "omegal.~a" i) ) 
(intern (format nil "omega2 . ~a" i) ) 
(intern (format nil "omega3.~a" i) ) 
(intern (format nil "omega4.~a" i) ) 
(intern (format nil "omega5.~a" i) ) 
(intern (format nil "omega6.~a" i))) 

(intern (format nil "a~a" i)) 

(format nil "Omega (~a)" i) 

(list (intern (format nil "theta.~a" i) ) 
(intern (format nil "theta.~a" (- i 1)))) 

) 

(send (symbol-value (intern (format nil "omega. "a" i))) :location 
(list (* (+ i 0.5) (/ 1 (+ plotnum 2))) 0.2) 

) 

; ; create the nu node 

(create-moral-node (intern (format nil "nu.~a" i)) 

(list (intern (format nil "nul.~a" i)) 
(intern (format nil "nu2.~a" i) ) 
(intern (format nil "nu3.~a" i) ) 
(intern (format nil "nu4 . ~a" i)) 
(intern (format nil "nu5.~a" i)) 
(intern (format nil "nu6.~a" i) ) ) 

(intern (format nil "b~a" i)) 

(format nil "Nu(~a>" i) 

(list (intern (format nil "theta.~a" i) ) 
(intern (format nil "x.~a" i))) 
) 

(send (symbol-value (intern (format nil "nu. "a" i))) :location 
(list (* (+ i 0.5) (/ 1 (+ plotnum 2))) 0.6) 

) 

; ; create the x node 

(create-moral-node (intern (format nil "x.~a" i) ) 



(list 


( intern 


(format 


nil 


"xl. "a 


( intern 


( format 


nil 


"x2 . 


"a" 


i) ) 


(intern 


( format 


nil 


"x3. 


"a" 


i) ) 


( intern 


( format 


nil 


"x4 . 


"a" 


i) ) 


( intern 


( format 


nil 


"x5 . 


"a" 


i) ) 


( intern 


( format 


nil 


"x6 . 


"a" 


i) ) ) 


(intern (format nil "b 


"a" 


i) ) 


(format nil "X(~a) 


" i) 






(list 


( intern 


(format 


nil 


"theta 


(intern 


(format 


nil 


"nu . 


"a" 


i) ) ) 



) 

(send (symbol-value (intern (format nil "x.~a" i) ) ) :location 
(list (* (+ i 1) (/ 1 (+ plotnum 2))) 0.8) 

) 
) 

The plots can now be created in the usual way. 

( create-moral-plot 'myplot) 
(create-global -moral -plot ' myplot 2 ) 
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3.6 Adjustments 



The first 6 weeks of observations will be added into the model. 



; ; Sequentially introduce the data 

(dolist (i (iseq 1 6)) 

(format t ""SData for week "a" i) 



(send (symbol-value 
(list 



(intern (format nil "b~a" i))) :observe 



( intern 


( format 


nil 


"xl . 


"a" 


i) ) 


( intern 


( format 


nil 


"x2 . 


"a" 


i) ) 


( intern 


( format 


nil 


"x3. 


"a" 


i) ) 


( intern 


(format 


nil 


"x4 . 


"a" 


i) ) 


( intern 


( format 


nil 


"x5 . 


"a" 


i) ) 


( intern 


( format 


nil 


"x6 . 


"a" 


i) ) 



(select (column-list data) (-11)) 
) 

(send myplot : record) 

(send (symbol-value (intern (format nil "b~a" 



i) ) ) : absorb) 



The resulting plot windows give a good impression of the adjustment process, and the way 
information flows forward and backwards through time in such models. 
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_ ++•#••• 

0mega(1)0mega(2)0mega(3)0mega(4)0mega(5)0mega(6)0mega(7)0mega(8) 

MMM \/\ 









Theta(O) Theta(1) Theta(2) Theta(3) Theta(4) Theta(5) Theta(6) Theta(7) Theta(8) 

/ 



/////// 






NuH) Nu(2) Nu(3) Nu(4) Nu(5) 





Nu(6) Nu?7) Nu?8) 



\ \ \ \ \ \ \ \ 






X(3) X(4) X(S) X(6) X(7) X(B) 
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Omega(1) Omega(2) Omega(3) Omega(4) Omega(5) Omega(6) Omega(7) Omega(8) 



Theta(O) Theta(1) Theta(2) Theta(3) Theta(4) Theta(5) Theta(6) Theta(7) Theta(8) 




Nu(1) Nu(2) Nu(3) Nu(4) Nu(5) Nu(6) Nu(7) Nu(8) 



* * * 



X(1) X(2) X(3) X(4) X(5) X(6) X(7) X(8) 



The file included as part of the distribution then goes on to extract the adjusted expectations of 
the residuals, and plot them on a line graph. See the example source file for more details. 



4 Important note/Disclaimer 

It is important to note that BAYES-LIN is a rapidly developing prototype system, and does contain 
many bugs. You should not rely on BAYES-LIN producing correct output, and should verify parts 
of calculations as far as possible, using alternative software, such as [B/D] (see [JTQ|]). The author 
accepts no liability whatsoever regarding the use of BAYES-LIN, errors or losses arising from the 
use of BAYES-LIN, etc. Feel free to check the source code, correct it, and email the corrections to 
the author. 
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5 Command reference 



This section lists all global variables, functions, object prototypes and methods defined by the 
BAYES-LIN system. On-line help is available. To obtain help on a global function, eg. create- 
moral -node use the expression (help ' create-moral-node) To obtain help for a method, 
eg. the : observe method of the moral-node-proto object, use the expression 

(send moral-node-proto :help :observe). 



5.1 Global variables 



Variable 


Description 


*tree-nodes* 
*moral-nodes* 


A list of symbols representing instances of tree-node objects cre- 
ated using the create-tree-node global function. 
A list of symbols representing instances of moral-node objects 
created using the create-moral-node global function. 



5.2 Global functions 



Function 


Description 


create-tree-node 
create-moral-node 
create-moral-plot 

ginv 


A function to create and initialise clique-tree objects. 
A function to create and initialise moral graph objects. 
A function to create a graphics window for illustrating and diag- 
nosing Bayes linear adjustments. 

Function to return the Moore-Penrose generalised inverse of a real 
square symmetric matrix. 



5.3 Object prototypes 



Object prototype 


Description 


tree-node-proto 
moral-node-proto 
moral -pi ot-proto 

gl obal -moral -pi ot-proto 


The prototype for objects representing clique- 
tree nodes. 

The prototype for objects representing moral 
graph nodes. 

The prototype for a graphics window object for 
the displaying of information relating to current 
Bayes linear adjustments. 
A plot to summarise the partial adjustments 
shown on the moral-plot-proto plots. 
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5.3.1 tree-node-proto slots and methods 



Slot 


Description 


name 
variables 
neighbours 

variance 
expectation 
var-d-inv 
cov-d-self 
obs-vars 
obs-d-ed 
location 


The name of the object. 

List of random variables associated with this tree-node object. 

List of neighbouring junction tree nodes. 

variance matrix associated wiin ine vanauie nsi. 

Expectation vector associated with the variable list. 

Var(D)- 1 

Cov (D, •) 

D 

d-E(D) 

On a (0, 1) x (0, 1) scale for plotting. 



Method 


Description 


: absorb 


Absorb information from last : observe ready for next observe. 


: aex 


Adjusted expectation. 


: avar 


Adjusted variance. 


: cov 


Current covariance. 


: ex 


Current expectation. 


: info 


Prints some information relating to the object. 


: location 


Accessor method. 


: observe 


Method to introduce data to the tree. 


: positions 


Variable positions. 


: propagate 


Method used to propagate information around the tree. 


: remove-neighbour 


Remove a neighbour from the list. 


: resolution 


Partial resolution for the current adjustment. 


: rvar 


Resolved variance matrix. 


: size-ratio 


Partial size-ratio for the current adjustment. 


: transform 


Partial resolution transform for the current adjustment. 


: var 


Current variance matrix. 



5.3.2 moral-node-proto slots and methods 



Slot 


Description 


name 
variables 
tree-node 

print-name 
neighbours 
location 
var_b_inv 

ex_b 
resolutions 
size-ratios 
global-size-ratios 


Name of the object. 

List of random variables associated with this moral node object. 
Name of a tree-node which contains all of the variables in this 
moral node. 

A string for plotting purposes. 

A list of neighbouring moral graph nodes. 

On a (0, 1) x (0, 1) scale for plotting. 

Inverse of the a priori variance matrix for the variables repre- 
sented by this node. 

The prior expectation vector for this node. 
List of resolutions for the partial adjustments. 
List of partial size-ratios for the adjustments. 
List of global size-ratios for the adjustments. 
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Method 


Description 


: absorb 


Absorb info ready for next : observe. 


: aex 


Adjusted expectation vector. 


: avar 


Adjusted variance matrix. 


: bearing 


Bearing vector. 


: ex 


Current expectation vector. 


: info 


Prints some info about the object. 


: location 


Accessor method. 


: observe 


Introduce data into the graph. 


: remove-neighbour 


Remove a node from the neighbour list. 


: resolution 


Partial resolution wrt a priori structure. 


: rvar 


Partial resolved variance matrix. 


: size-ratio 


Partial size-ratio. 


: global-size-ratio 


Global size-ratio. 


: transform 


Partial resolution transform. 


: tree-node 


Accessor method. 


: var 


Current variance matrix. 



5.3.3 moral-plot -proto slots and methods 

This object inherits all slots and methods from graph-window-prot o, but also has the follow- 
ing. 



Slot 


Description 


nodes 
real-size 
radius 
diagnostics 
node-labels 
outlines 


List of nodes to be plotted. 

Window size. 

Node radius (scaled). 

Flag for diagnostic plotting. 

Flag for node label printing. 

Flag for node outline printing. 




Method 


Description 


: plot-arcs 
: plot-node 
: r-to-s 
: record 
: redraw 
: resize 
: s-to-r 
: sr-map 

: diagnostics 
: node-labels 
: outlines 


Draw the arcs associated with a given node. 

Draw the given node, and all its shadings. 

Take "real" (screen) coords to scaled coords. 

Record current adjustment information for inclusion in the plot. 

Guess! 

Recalculate scale parameters. 
Scaled to real coord transform. 

Function which maps (0,°°) to ( — 1, 1) monotonically, mapping 1 

to 0. This is the function used to transform size ratios for red and 

blue diagnostic shadings. 

Set and unset diagnostics plotting. 

Set and unset node label printing. 

Set and unset node label printing. 



5.3.4 global-moral-plot -proto slots and methods 

This object inherits all slots and methods from moral-plot-proto and has no others. 
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